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Abstract 

New measurements of the spin structure functions of the proton and deuteron gf(x,Q 2 ) and 
gf(x,Q 2 ) in the nucleon resonance region are compared with extrapolations of target-mass- 
corrected next-to-leading-order (NLO) QCD fits to higher energy data. Averaged over the entire 
resonance region (W < 2 GeV), the data and QCD fits are in good agreement in both magnitude 
and Q 2 dependence for Q 2 > 1.7 GeV 2 /c 2 . This "global" duality appears to result from cancel- 
lations among the prominent "local" resonance regions: in particular strong a^/ 2 contributions in 
the A(1232) region appear to be compensated by strong 0\i<i contributions in the resonance region 
centered on 1.5 GeV. These results are encouraging for the extension of NLO QCD fits to lower W 
and Q 2 than have been used previously. 

PACS numbers: 13.60.Hb, 25.30.Fj,24.30.Gd 
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The theoretical description of particle interactions has utilized quark-gluon degrees of 
freedom at high energies and hadronic degrees of freedom at low energies. With suitable 
averaging over resonant excitations, the two approaches have been found in several cases 
to be nearly equivalent, a phenomenon referred to as quark-hadron duality. These cases 
include e + e~ annihilation, semi-leptonic decays of heavy mesons, electron-pion scattering, 
semi-inclusive deep-inelastic scattering, and both spin-averaged and spin- dependent inclu- 
sive lepton-nucleus scattering the subject of the present investigation. Pragmatically, 
understanding the limitations and applicability of quark-hadron duality in this process is 
useful in order to define the kinematic region in which parton distribution functions (PDF) 
can be reliably extracted. 

In lepton-nucleon scattering, the low and high energy regimes have traditionally been 
separated using W, the invariant mass of the hadronic final state, and Q 2 , the four mo- 
mentum transfer squared. A region of prominent nucleon resonances is observed for W < 2 
GeV and Q 2 < 10 GeV 2 /c 2 , while for higher W or Q 2 there is no longer any obvious res- 
onance structure. Historically, quark-hadron duality was first observed in 1970 by Bloom 
and Gilman |2j in the spin-averaged lepton-nucleon process. They noted that the inclusive 
structure function, F 2 (W,Q 2 ) averages smoothly at low W and Q 2 to the scaling function 
Fz(W, Q 2 ) measured at high energy, using an empirical scaling variable in place of the original 
Bjorken x scaling variable. Subsequently, Georgi and Politzer found that quark-hadron 
duality is exhibited down to Q 2 ~ 1 GeV 2 /c 2 using the Nachtmann [4] scaling variable 



£ = 2x/(l + a/1 + 4M 2 x 2 / Q 2 ) which approximates the purely kinematic higher twist cor- 
rections arising from the non-zero nucleon mass M. More recently, explicit target-mass 
(TM) corrections have been derived in the framework of QCD for both unpolarized and 
polarized structure functions that obviate the need for an approximate scaling variable. 

To explain quark-hadron duality theoretically, de Rujula, Georgi, and Politzer [f| em- 
ployed a perturbative operator product expansion of QCD structure function moments. In 
this framework, quark-hadron duality implies a small net effect from higher twist contri- 
butions, once the kinematic target-mass contribution is taken into account. In a simple 
QCD picture, the additional higher twist contributions (which are proportional to powers of 
1 / a/Q 2 ) are due to quark-quark and quark-gluon correlations. Close and Isgur provided 
an interesting explanation in the constituent quark model in terms of cancellations from 
resonance contributions with opposite parity. A recent theoretical QCD study 8] of both 
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polarized and unpolarized structure functions incorporates many of these concepts, with the 
addition of a careful treatment of so-called high-x resummation corrections. The authors 
concluded that higher-twist corrections are suppressed more for the unpolarized structure 
function F 2 than for the proton polarized structure function g\, where a sizable negative 
contribution is observed. A comprehensive review of quark-hadron duality from both the 
experimental and theoretical perspective was also published recently [l|. 

Unpolarized structure function data exhibit excitation-energy-averaged scaling not only 
averaged over the entire resonance region (M < W < 2 GeV), referred to as "global duality", 
but also in each of several restricted regions in W, corresponding to the three prominent 
resonance regions centered on W = 1.23, 1.5 and 1.7 GeV, a phenomenon referred to as 
"local duality" . This was demonstrated experimentally using high-statistical-accuracy data 
from Jefferson Lab [lOj, and interpreted theoretically by Carlson and Mukhopadhyay 
using the expected pQCD Q 2 -dependence of nucleon transition form factors. 
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The new data presented in this report augment previously available results for g\ from 
SLAC DESY [3| and JLab [3, [lj| with higher statistical precision and an expanded 

range of Q 2 . This allows us to experimentally examine local duality for g\ much more 
accurately than was previously possible. The addition of a considerable body of deuteron 
gf data allows the first examination of the isospin dependence of global duality in g\. 

When testing duality there is an intrinsic uncertainty as to which DIS curves to use for 
comparison with the averaged resonance region data. In this paper, we choose the average 
of two representative Next-to-Leading Order (NLO) QCD fits 13, to polarized structure 
function data above the resonance region. The NLO evolution is considered to be reasonably 
reliable down to Q 2 values of order 1 GeV 2 /c 2 . We choose to use fits with NLO evolution, 
rather than LO or purely empirical fits to data, to give the best possible estimate of the 
Q 2 -dependence of g±. The high-energy data used in the NLO QCD fits have relatively large 
errors compared to those for unpolarized structure functions, particularly at high values 
of x that tend to correspond to our resonance region data. Since precise error bands in 
our kinematic region are not available, we ascribe a very approximate relative error of 10% 
(20%) to the g\ (gf) DIS fits, independent of x. The error on the average deuteron DIS fit 
is larger due to the much larger relative contribution of negatively polarized quarks in the 
neutron compared to the proton. We take kinematic target-mass corrections into account 
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using the prescription of Bliimlein and Tkabladze 



(x + £)j f 1 du 



[ -9? CD (u,Q> 



^(l + 7) 2 J* 
7(2-7) f'du rdVgQCD^ 



2(1 + 7)5/2 u J u V 

where 7 = 4M 2 x 2 /Q 2 . This prescription is not unique, and in particular has the drawback 
of resulting in non-zero values of g% at x — 1. An approach that avoids this problem has been 
worked out for F2 [3], but is not yet available for g\. The calculation of high-x resummation 
corrections is theoretically more complicated [^| and technically more challenging than target 
mass corrections. Rather than attempting these calculations ourselves, we simply note that 
Ref. H finds enhancements of order 10% to 20% for the proton averaged over the full 
resonance region, roughly independent of Q 2 for 0.5 < Q 2 < 5 GeV 2 /c 2 . This corrections 
could well be different for the deuteron and for individual "local" regions in W. 

The ana^is is based on recent* pnbbshed data Q from Jefferson Lab. Very briefly, to 
this experiment the CEBAF Large Acceptance Spectrometer 2l( in Jefferson Lab's Hall B 
was used to measure spin asymmetries in the scattering of longitudinally polarized electrons 
from longitudinally polarized protons and deuterons. The data were collected in 2001 using 
incident energies of 1.6 GeV and 5.7 GeV. Beam currents ranged from 1 to 5 nA, and the 



beam polarization averaged 70%. The detector package [2l( allowed clean identification 
of electrons scattered at polar ang les between 8 and 45 degrees. Ammonia, polarized via 
Dynamic Nuclear Polarization [22j, was used to provide polarized protons and deuterons, 
using the 15 NH3 and 15 ND3 isotopes, respectively. The average target polarization was about 
75% for the proton and about 25% for the deuteron. The data were divided into 40 bins in 
Q 2 , equally spaced on a logarithmic scale between 0.01 and 10 GeV 2 /c 2 . 

Values of gi(x, Q 2 ) were determined from the ratios of gx/F\ presented in j^l using recent 



fits to proton 



2311 and deuteron 



24j data to evaluate the unpolarized structure function 
Fi(x, Q 2 ). The resulting values of g±(x, Q 2 ) for both the proton and the deuteron are plotted 
(scaled by x) as a function of x for four representative Q 2 bins in Figure ^ The three arrows 
on each panel correspond to the three prominent resonance regions at W = 1.7, 1.5, and 
1.23 GeV, from left to right. We compare the data to the extrapolations of DIS fits (as 
described above), represented by the hatched bands. 



It can be seen in Fig. 1 that the data in fixed Q 2 bins indeed exhibit oscillations in x 
compared to the smooth behavior of the DIS curves. In addition, the averaged data become 
increasingly commensurate with the models with increasing Q 2 , as expected if quark-hadron 
duality is valid for g\. In closer detail, one can also observe that the experimental data for 
both the proton and the deuteron lie below the curves in the A (1232) region, and above 
in the W = 1.5 GeV [S u ( 1535) /D 13 ( 1520) /P u ( 1440)] region. This is not surprising at low 
to moderate Q 2 , where resonant contributions dominate over non-resonant contributions. 
Recall that g\ is proportional to 01/2 — 03/2- The N — > A(1232) is known to be dominated 
by Ml strength [25j over the Q 2 range of the present study, which results in a virtual photon 
cross section cr 3 / 2 about three times larger than 01/2, corresponding to negative values of 
Qi. In contrast, in the DIS limit of incoherent scattering from massless quarks, g\ must 
be positive due to helicity conservation. In the second resonance region, the «Sii(1535) 
and Pn(1440) transitions can only contribute to 01/2, and recent studies show that 
the N — > Di3(1520) changes from dominantly 03/2 at low Q 2 to dominantly 01/2 above 1 
GeV 2 /c 2 . The three resonances together therefore are expected to result in large positive 
values of g\ above 1 GeV 2 /c 2 (potentially larger than the DIS limit). 

To clarify these observations with respect to both local and global duality, we have 
averaged over x both data and models for gi over a Q 2 -dependent interval corresponding to 
four specific regions in W. The x-averaged values of gi for the entire resonance region (scaled 
by Q 2 ) are plotted as a function of Q 2 in Fig. [2] for both targets. The proton averages for 
four smaller regions in W are plotted in Fig. 01 Specifically, the averages were determined 
as 



X 



<9i(Q 2 )>-~ 

Xfo X\ 

where Xi and Xh correspond respectively to the maximum and minimum values of W in the 
interval considered, at the given value of Q [using the definition x~ l = 1 + (W 2 - M 2 )/Q 2 }. 
The TM-corrected NLO PDF parameterizations shown in Fig. 1 were averaged in the same 
way as for the experimental data. 

The averages displayed in Figs. 2 test "global" duality by averaging g\ over x for the 
entire region from pion threshold to W = 2 GeV. The data for both targets exhibit a power- 
law-type deviation from the DIS curves at low Q 2 , but essentially agree with them above 
Q 2 = 1.7 GeV 2 /c 2 , within the systematic errors of the data and models. This is somewhat 
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higher in Q 2 than for the unpolarized F2 structure, supporting the conclusions of Ref . |8| . 

Turning to the examination of "local" duality for the proton, it can be seen in the upper 
left panel of Fig. 3 that in the "first" resonance region, dominated by the A(1232) resonance, 
the data have the opposite sign of the extrapolations of DIS models at low Q 2 . This is allowed 
by the spin-3/2 nature of the A(1232), and is expected due to the dominance of the Ml 
transition strength [25[, as discussed above. What is interesting is that, while the data 
change sign at Q 2 near 1 GeV 2 /c 2 , the averaged values are significantly below the models 
even to the highest Q 2 of the present experiment, in spite of the fact that the N — > A(1232) 
transition form factor (FF) decreases more rap idly with Q 2 than, for example, the elastic 



FF or the N 5 U (1535) transition FF |2j,|27| (a phenomena sometimes referred to as the 
"disappearing A"). It is evident that the A(1232) has not yet completely disappeared at 
Q2 = 5 GeV 2 /c 2 . 

In the second resonance region, two of the three known resonances [Pu(1440) and 
Sii(1535)] contribute only to am, while the third [_Di3(1520)] contributes more to cryz 
than to (T3/2 above 1 GeV 2 /c 2 j2q. Therefore, it isn't surprising that the proton data lie 
significantly above the DIS extrapolations in this narrow region of W . 

In the "third" resonance region centered on 1.7 GeV, the F 15 (1680) is dominant at low Q 2 , 
but above about 1 GeV 2 /c 2 other resonances are also important [26j, such as the Sii(1650), 
631(1620), and D 33 (1700). The F 15 contributes mainly to a 3 / 2 at low Q 2 (i.e. negative gi), 
but switches to 01/2 dominance at higher Q 2 [26] . The average over all of these resonances 
plus non-resonant background produces very good agreement between data and DIS models 
in this region, as might be expected from the parity-averaging arguments of Close and 
Isgur 

For completeness, we have also studied a fourth region centered on 1.9 GeV (for which 
there are a large number of poorly established resonances, difficult to distinguish from non- 
resonant contributions). In this case, the data lie slightly below the DIS models, although 
the significance is marginal when systematic errors are taken into account. It appears that 
much of the good agreement between data and models observed in the entire resonance 
region comes about from pairing the "first" and "second" resonance regions together, with 
further improvement from including the "fourth" region. This lends further support to the 
Close-Isgur model. 

Again following Close and Isgur 7], one might expect DIS and resonance region data to 
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converge at lower values of Q 2 if the ground state elastic contribution is also included in 
the global duality averaging. The open circles in Fig. 2 include the elastic (quasi-elastic) 
contributions to the g\ averages for the proton (deuteron), given by G e(Q 2 )[G e(Q 2 ) + 
tGm(Q 2 )]/2(1 + r)(xh — Xi), where r = Q 2 /4M 2 . To evaluate the nucleon electric and 
magnetic form factors G r(Q 2 ) and Gm(Q 2 ), we used a slightly modified version of the 
parameterization of Ref. [28(. For both the proton and the deuteron, the Q 2 - dependence 
with the elastic contribution more closely resembles the Q 2 -dependence of the models, down 
to values of Q 2 as low as 0.7 GeV 2 /c 2 , which is already pushing below the expected region 
of validity of 1 GeV 2 /c 2 for the PDF fits. However, the magnitude of the data is of order 
10% to 20% higher than the models. As mentioned above, high-x resummation corrections 
may well account for this difference. The result of pairing the A(1232) resonance (the lowest 
spin-3/2 ground state) with the elastic contribution (the lowest spin-1/2 ground state) is 
illustrated in the upper left panel of Fig. 3 for the proton. The elastic contribution, in 
the way we have treated it, over-compensates, resulting in power-law deviations at low Q 2 
that lie well above the data, rather than well below. It appears that including the elastic 
contribution with the entire resonance region works much better than pairing it with the 
single A(1232) resonance. 

In summary, we have used data for both the proton and deuteron to examine both "local" 



and "global" quark-hadron duality in g\. As was determined in previous studies [l|, Isj], g\ in 
the resonance region oscillates around extrapolations of NLO PDF fits to higher energy data, 
especially when target-mass corrections are taken into account. Averaged over the traditional 
resonance region (W < 2 GeV), the data and fits agree within errors for Q 2 > 1.7 GeVVc 2 , 
a slightly higher value than observed for the spin-averaged structure function F 2 ^, llo| . 
Including the elastic contribution may extend the region of agreement to Q 2 = 0.7 GeV 2 /c 2 , 
after consideration of the uncertainties due to high-x resummation. A similar effect was 
found in the unpolarized case . We find similar results for the previously unexamined 

gf structure function as for g\, indicating no large effects from different isospin projections. 
In terms of "local" duality, we find that in the A(1232) region, the proton data lie below the 
PDF fits, even for Q 2 values as large as 5 GeV 2 /c 2 , while the region centered on 1.5 GeV 
lies above the PDF fits for all Q 2 values studied. It appears that global duality is largely 
realized by summing over the four lowest-mass resonances. Since the Q 2 - dependence for the 
PDF fits and the data are remarkably similar above Q 2 = 1.7 GeV 2 /c 2 , we conclude that 
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from the practical point of view, it is not unreasonable, with suitable averaging over W, to 
include data with W > 1.58 GeV and Q 2 > 1.7 GeV 2 /c 2 in future global NLO PDF fits, 
as long as the effects of TM and high-x resummation effects are taken into account. In the 
near future, high statistical accuracy data from the present experiment with 2.5 and 4.2 
GeV electron energies, and additional data at 5.7 GeV, presently under analysis, will allow 
more precise studies, particularly for the deuteron. 
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FIG. 1: Present data for proton g±(x,Q 2 ) (left panels) and deuteron gf(x,Q 2 ) (denned to be 
per nucleon, right panels) at four representative values of Q 2 . The errors include statistical and 
systematic contributions added in quadrature. The three arrows on each plot indicate the central 
kinematic position of the three prominent resonance regions at W = 1.7, 1.5, and 1.23 GeV from 
left to right. The hatched band represents the range ofoi predicted by modern NLO Parton 
Distribution Function (PDF) fits (GRSV 3 

and AAC [13]) to high energy data, evolved to the 



Q 2 of our data and corrected for target-mass as described in the text. 
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FIG. 2: The Q 2 -dependence of Q 2 gi(x,Q 2 ), averaged over a region in x corresponding to 1.08 < 
W < 2 GeV (solid circles) for: a) proton; b) deuteron. The inner error bars reflect only statistical 
contributions, while the outer error bars include statistical and systematic components added in 
quadrature. The open circles represent our data after adding the contribution from ep elastic (ed 
quasielastic) scattering at x = 1 for the proton (deuteron). For clarity these results are slightly 
displaced in Q 2 , and the error bars include only statistical contributions. The hatched bands 
represent the range of the averages calculated from extrapolated NLO DIS fits, as in Fig. 1 (see 
text for details). 
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FIG. 3: The Q 2 -dependence of Q 2 g±(x, Q 2 ) for the proton, averaged over various regions of x. At 
each Q 2 , the x-range over which gi was averaged is determined by the corresponding range in W 
as indicated in each panel (see text). Symbols and curves as in Fig. [3 
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